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A -5 M ffi SUPER-EARTH ORBITING GJ 436?: THE POWER OF NEAR-GRAZING TRANSITS 
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ABSTRACT 

Most of the presently identified exoplanets have masses similar to that of Jupiter and therefore are 
assumed to be gaseous objects. With the ever-increasing interest in discovering lower-mass planets, 
several of the so-called super-Earths (1 M®< M <10 M®), which are predicted to be rocky, have 
already been found. Here we report the possible discovery of a planet around the M-type star GJ 436 
with a minimum mass of 4.7±0.6 M® and a true mass of ~5 M® , which would make it the least massive 
planet around a main-sequence star found to date. The planet is identified from its perturbations on an 
inner Neptune-mass transiting planet (GJ 436b), by pumping eccentricity and producing variations in 
the orbital inclination. Analysis of published radial velocity measurements indeed reveals a significant 
signal corresponding to an orbital period that is very close to the 2:1 mean motion resonance with the 
inner planet. The near-grazing nature of the transit makes it extremely sensitive to small changes in 
the inclination. 

Subject headings: Stars: planetary systems — planetary systems: formation — stars: individual 
(GJ 436) 



1. INTRODUCTION 

Hundreds of exoplanets have been discovered over the 
past decade using a variety of techniques, notably pre- 
cision radial velocities, most of them being Jupiter-like. 
While such planets are extremely useful to understand 
the morphology and evolution of planetary systems and 
star-planet connections, there is an obvious interest in 
eventually identifying an Earth-like object. Model calcu- 
lations indicate that planets with masses in t he interval 
1 Mffl< M <10 M m are of terrestrial type (e.g. Jlda fc Linl 
120041: 1 Valencia et all2007f) . Such planets have often been 
dubbed "super-Earths". In this quest for smaller plan- 
ets there have already been some successful detections of 
planets that fall in the super-Earth domain (|Rivera et al.l 
120051: iBeaulieu et~all 120061: lUdrv et alj 120071), arid ev en 
their habitability may be evaluated (jSelsis et al.ll2007j ). 

Besides the direct detection from radial velocities or 
transits, methods for the discovery of low-mass planets 
as perturbers to other planets have been suggested (al- 
though without any successful detections yet). These 
are mostly based on monitoring variations in the timing 
of the transit over a relatively long time scale so that 
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In this Letter we make use of yet another method that 
has allowed the possible detection of a super-Earth per- 
turbing the transiting planet GJ 436b. In this case, it 
is not perturbations to the transit mid-time but to the 
overall orbital elements of the inner transiting object that 
cause variations on the transit shape and depth. We show 
that, for near-grazing transits, their duration is strongly 
sensitive to small changes in the orbital inclination. 

2. AN UNUSUAL HOT-NEPTUNE AROUND GJ 436 
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The M2.5-dwarf GJ 436 was discovere d to host a 
Neptu ne-mass planet in a 2.6-d orbit by iButler et al.l 
(2004). Two properties made this object especially in- 
teresting, namely its relatively small mass and a surpris- 
ing non-zero eccentricity of about 0.15. Such value of 
the eccentricity was re cently confirmed by the analysis of 
IManess et all (|2007t l. IButler et all (|2004h also obtained 
high-precision photometry to investigate the presence of 
transits but ruled out the possibility of a transit with 
a depth greater than 0.4%. However, a surprise came 
with the actual detection of transits f rom GJ 436b with a 
depth of 0.7% bv lGillon et~all (|2007bD . thus becoming, by 
far, the smallest transiting planet yet detected. A series 
of studies, mostly using Spitzer Space Telescope (SST), 
have greatly c ontributed to establ i shing the properties 
of the planet (IDeming et all 120071 : iGillon et all |2007at 
iDemorv et al.l feOO?: Torres 2007) and also to strengthen 
the case for an eccentric orbit by observing the occulta- 
tion event at orbital phase 0.59. 

The origin of the high eccentricity o f GJ 436b 
was investigated in detail by IManess et al.l (|2007l ) and 
IDeming et al.l (|2007l ) . Both studies conclude that the cir- 
cularization timescale (~10 8 yr) is significantly smaller 
than the old age of the system (>6T0 9 yr) when assum- 
ing reaso nable values fo r the p lanet's tidal dissipation pa- 
rameter. IManess et al.l (|2007l ) also pointed out the pres- 
ence of a long-term trend with a value of 1.3 m s _1 per 
year on the systemic radial velocity of GJ 436. Thus, the 
authors investigated the possibility that the eccentricity 
and the long-term velocity trend could be explained from 
the perturbation exerted by an object in a wider orbit 
without reaching conclusive results. 

We propose an alternative possibility to explain the 
eccentricity of GJ 436b, namely the perturbation from 
a relatively small planet in a close orbit. We show be- 
low that the effects of the perturber on the inner planet 
can excite the eccentricity up to the observed value. But 
G J 436b has yet another characteristic that makes it dif- 
ferent to other transiting planets and this is the near- 
grazing nature of its transit. The impact parameter of 
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the transit was found to be about 0.85, which implies an 
orbital inclination of 86?3. If the inclination happened 
to be just 85?3 the planet would not cross the disk of the 
star. Studies mostly focused on triple star systems have 
pointed out that perturbations may not only change the 
eccentricity but also other orbital properties of the in- 
ner orbit. When the perturber resides in a non-coplanar 
orbit this gives rise to a modulation in th e inclination 
of the inner object (e.g., ISoderhielmlll975l ). As an ex- 
ample, the triple system scenario has been advocated to 
explain the cessation of eclipses in th e binary SS Lac 
(|Torres fc Stefanikl 120001 iTorred [200lh. In the contex t 
of exoplanets, [ Schncidci] (|1994l ). iMiralda-Escudl (|2002f ). 
lLaughlin et ali i 20051 ) considered the possibility that pre- 
cession induced by a perturbing planet could lead to ob- 
servable changes in the duration (and existence) of the 
transit of an inner giant planet. 

GJ 436b makes an ideal system to find evidence for a 
perturbing small planet, because of the telltale non-zero 
eccentricity, but also to put severe constraints on the 
properties of the perturber owing to the extreme sensi- 
tivity of the current configuration to small changes in the 
orbital inclination angle. 

3. A SECOND PLANET AROUND GJ 436? 
3.1. Dynamical study 

A possible explanation to the apparently contradict- 
ing results concerning the detection of transits is that 
the orbital inclination has indeed changed during the 
3.3-year interval between the different photometric ob- 
servations. Calculations show that an orbital inclina- 
tion <86° w o uld h ave made the transit undetectable to 
iButler et al.l (|2004h 's photometric measurements. From 
these considerations a small variation of the inclination 
ang le at a rate of roug hly ~0?1 yr" 1 coul d make both 
the IButler et al.l (|2004l ) non detection and iGillon et al.l 
(|2007bl )'s discovery of transits compatible. Note that 
this is only a possible scenario since the photometry of 
IButler et al.l has relatively sparse phase coverage. In ret- 
rospect, from the currently known durati on and prop- 
erties of the transit, 3 measurements from IButler et ahl 
(2004) should have betrayed the presence of the planet, 
although with low significance. 

Assuming this hypothesis, it is reasonable to explore 
the possibility of a perturber that could be responsible 
for both the relatively large eccentricity and the incli- 
nation change, while remaining undetected by the ra- 
dial velocity measurements. From this, one can assume 
that the semi-amplitude of the perturber should be be- 
low about 4 m s _1 , implying that its mass should satisfy 
the following inequality: M p < 30 a 1 / 2 M®, with a being 
the orbital semi-major axis in AU. On the other hand, 
because of the inclination change, the variation of the 
transit duration with time is of about 10~ 5 . Following 
iMiralda-Escudel (|2002f ) we find that the perturber mass 
should satisfy: M p > 3 • 10 4 a 3 M®. From both inequal- 
ities, there is an allowed range of perturber masses and 
semi-major axes. 

For more accurate estimates we carried out direct in- 
tegration s of the equation s of motion using the Mercury 
package (jChambersI Il999f ) . We started with an inner 
planet in a circular orbit and with the currently observed 
semi-major axis. Then, we considered different combina- 
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Fig. 1. — Allowed region for a perturbing planet to GJ 436b 
given the observational constraints. The dashed and dotted lines 
correspond to the radial velocity detection limits for inclinations 
of 90° and 45°, respectively, and the solid line represents the limit 
from an approximate calculation for a perturber (see text). The 
shaded regions illustrate the allowed range of masses and semi- 
major axes for the perturbing object. The circles correspond to 
a grid search to identify cases in which a planetary companion to 
GJ 436b could be able to induce the observed eccentricity and rate 
of orbital inclination change. Positive cases are marked with filled 
circles while empty circles indicate configurations with negative 
results. The square marks the position of the planet detected from 
the radial velocity analysis. 

tions of mass (from 1 to 14 M®), semi- major axis (from 
0.04 to 0.1 AU), eccentricity (from 0.05 to 0.3) and in- 
clination (from 85° to 45°) for the perturber. The inte- 
grations were performed for a time interval of 10 5 yr to 
guarantee the stability of the planetary systems. In Fig. 
[T] we provide an illustration of the region in the mass 
vs. semi-major axis plane where perturbers can meet all 
constraints, i.e., sufficient eccentricity of the inner planet 
and minimum inclination variation rate. 

We further explored semi-major axis values at mean- 
motion resonances (MMRs). Location in a MMR can 
be a stabilizing factor and also perturbations ca n reach 
their maximum efficiency (e.g., Ag ol et all 120051 ). Inte- 
grations for semi-major axes corresponding to the fol- 
lowing MMRs were carried out: 3:2, 5:3, 2:1, 3:1, and 
4:1. In all cases, the presence of the planet in a MMR 
increased the stability and, further, perturbing planets 
with smaller masses were able to induce the observed 
eccentricity and orbital inclination change to the inner 
planet. For the strongest 2:1 resonance we found a lower 
limit to the perturbing planet mass of only 1 M® at an ex- 
treme eccentricity and relative inclination. For the case 
of a perturbing planet with 3-7 M®, eccentricity values 
of 0.15-0.20 and initial inclination differences of only 5- 
15° were sufficient to explain the observed eccentricity 
and rate of inclination change of the inner planet. 

In the analysis we neglected tidal dissipation since we 
focus on the current snapshot of the orbital configuration 
of the system. The planets must be undergoing signifi- 
cant tidal dissipation because of the non-zero eccentric- 
ity. Our calculation of the dynamical evolution of the 
system is valid to first order because the tidal evolution 
timescale (^10 8 yr) is long compared with the timescale 
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of the orbital perturbations (~10 2-3 yr). Other effects 
have been neglected at this stage, which include preces- 
sion caused by the quadrupole moment of the star and 
by General Relativity. A more detailed dynamical study 
of the system is left for a future paper. 

3.2. Re-analysis of the radial velocity data 

To place more stringent constraints on the perturb- 
ing planet we studied the available radial velocity ob- 
servations, coveri ng a time lapse of 7 ye ars. The or- 
bital s olutions of iButler et al.1 (|2004f ) and iManess et al.1 
(2007) yield radial velocity semi-amplitudes of about 
18 m s _1 for GJ 436b, in a solution with a rms resid- 
ual of 4.8 m s _1 . We further analyzed the data by pre- 
whitening of the frequencies from the reported Neptune- 
mass planet. The periodogram revealed a relatively 
strong peak corresponding to a period of about 5.2 d. 
To evaluate the significance of the detection, we cal- 
culated the false-al arm probabil it y foll owing the Monte 
Carlo procedure in IButler et all ((20041 ). From 1000 re- 
alizations we conclude that the false-alarm probablility 
of the observed peak is ~20%, considering that only ob- 
jects with periods between ~4 and ~8 days could be 
responsible for the observed perturbations in the orbital 
eccentricity and inclination of th e inner planet. In a re- 
cent paper, iDemorv et al.1 (|2007f ) report 23 spectroscopic 
measurements using the HARPS sp ectrograph of GJ 436 . 
Combining these with the data from lManess et aTl (|2007l ) 
could provide the needed stronger proof of the signal 
studied here. 

Considering the evidence from the dynamic integra- 
tions, the peak in the periodogram is sufficiently signif- 
icant to merit analysis. Starting from the 5.2-d period, 
a simultaneous fit to the orbits of two planets yielded 
the results in Fig. [2] and Table [TJ In the fit we fixed 
the period of the transiting planet to 2.643913 d, which 
provides a good match to the recent ground-based transit 
timings, and the eccentricity and the argument of perias- 
tron to 0.15 and 343°, respectively, for co nsistency with 
the observation of the oc cultation event (jDeming et al.l 
120071: IDemorv et aHl2007[ ). Certainly, the fit to the ra- 
dial velocities should be consistent with our dynamical 
analysis. This is somewhat hampered by two uncon- 
strained parameters, namely the relative inclination of 
the planets and also the true orbital orientation (radial 
velocities are only sensitive to the argument of the peri- 
astron). We evaluated a range of possible rates of vari- 
ation in the eccentricities, arguments of periastron and 
orbital periods of the two planets and found them to 
have negligible effects on the fits given the uncertainty 
of the radial velocity measurements and the time base- 
line. However, we considered a fit allowing for a linear 
change in these three elements for the transiting planet 
and yielded e = 0.03 ± 0.02 yr" 1 , w = -1.8 ± 1.5° yr" 1 
and P — — 5.7±2.0 s yr -1 , all with low significance. The 
orbital elements in Table[T]should be regarded as osculat- 
ing elements at the mean epoch of the radial velocities. 
Note that the period given is not sidereal but anomalistic 
(i.e., apparent). From the latter (P), the sidereal period 
(P s ) can be computed as P s = P(l — wP/2ir). 

The two-planet fit reduces significantly the rms resid- 
uals of the radial velocities to 3.4 m s . We fixed the 
orbital eccentricity of the second planet because of insta- 
bilities in the solution. The adopted eccentricity of 0.2 is 




Orbital phase 

Fig. 2. — Two-planet radi al velocity fit to GJ 436. Radial ve- 
locity observations of GJ 436 (Mancss et al. 2007) were fitted with 
a model considering the orbital motions of two planets combined 
plus a long term radial velocity drift. The panels show the radial 
velocities (with la error bars) associated to each respective planet 
where the contribution from the other planet has been removed, 
together with the best orbital fit. 



TABLE 1 

TWO-PLANET FIT TO THE RADIAL VELOCITIES. 



Parameter 


GJ 436b 


GJ 436c 


P 1 (days) 
T pcri (HJD) 

e 

u (°) 

K (m s- 1 ) 
a (AU) 
A/sini (M®) 


2.643913 (fixed) 
2451551. 65±0.01 
0.15 (fixed) 
343 (fixed) 
18.4±0.4 
0.0287±0.0003 
23.2±0.5 


5.1859±0.0013 
2451553. 2±0.7 
0.2 (fixed) 
265±43 
3.0±0.4 
0.0450±0.0004 
4.7±0.6 


Radial velocity drift 
rms (m s _1 ) 


1.1±0.2 
3.36 
3.1 





Anomalistic period. 



compatible with our perturbation analysis and provides 
a good fit to the data. The planet's minimum mass is 
4.7±0.6 M®. From the perturbation analysis and plane- 
tary system formation arguments, the relative inclination 
of the two planets is likely to be below 15°, and therefore 
the real mass of the planet should be about 5 M ffi with 
12% uncertainty. While this could not be considered an 
extremely solid detection, its is significant and the planet 
has the correct properties to explain the orbital effects 
observed in GJ 436b. 

From the final system configuration, our dynamical cal- 
culations predict librations of the orbital elements of both 
planets over different timescales and with different ampli- 
tudes. Focusing on the inner planet, the inclination has a 
libration amplitude of ~5° with a period of ~110 yr. Its 
argument of periastron precesses with a period of ^500 yr 
but also has a libration amplitude of ~20° over a period 
of ~70 yr. For comparison, the General Relativistic pre- 
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Fig. 3. — Effect of a +0?1 change in the orbital inclination of 
GJ 436b. The upper panel shows t he best fit to the 8 fim transit 
data from SST IIGillon et al.ll2007aT) . The middle panel shows the 
residuals from the fit, which result in a x 2 value of 486 with 355 
degrees of freedom. The bottom panel illustrates the residuals of 
the fit when assuming an orbital inclination higher by +0?1, with 
a x 2 increasing to 578. 

cession period is ^15,000 yr. Finally, the eccentricity has 
a libration amplitude of ~0.1 with a period of ~70 yr. 

GJ 436c is the least massive planet known to orbit a 
main-sequence star and only the second bona-fide warm 
super-E arth together with t he 7.5±0.7 Mg planet around 
GJ 876 dl lvera since only minim um masses 

are a vailable for the planets around Gl 581 (|Udrv et al.1 
|2007| ). Interestingly, GJ 436c is found at an orbital pe- 
riod that is very close to the 2:1 MMR with GJ 436b, but 
not exactly so: P c /P b = 1.9614 ± 0.0005. This small but 
significant departure may be a consequence of the tidal 
evolution of the system. 

4. PROSPECTS FOR CONFIRMATION 

Because of the inclination change, the effects of the 
perturbing planet will become evident in high-precision 
transit photometry collected during 2008. We carried 
out tests to assess the capability to detect small changes 
in the inclination using real SST IRAC primary transit 
observations. We adopted a non-linear limb darkening 
law model with 4 coefficients and f ollowed the procedure 
described in lMandel fc Agoll (|2002D to fit the transit light 
curve. Then, we fitted the same dataset with models for 
an inclination larger by 0?1. Such expected inclination 
change will increase the transit duration in ~2 minutes 
making it detectable from SST at reasonably high confi- 
dence level. Fig. [3] shows the best fit to the 8 /im SST 
data and the residuals. 

5. CONCLUSIONS 

The proposed existence of GJ 436c is based on three 
independent pieces of evidence: a) An orbital inclina- 
tion change supported by the lack of transit detection 
in 2004 and the grazing transit observed in 2007; b) an 
inner planet with significant orbital eccentricity and a 
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tidal dissipation timescale much smaller than the age of 
the system; and c) a low-amplitude radial velocity sig- 
nal that is fully consistent with a planet in the 2:1 mean 
motion resonance with the inner planet. A key element 
of our study is a dynamical investigation that explains 
the observed effects a) and b), and predicts the existence 
of a perturbing planet with constraints on its mass and 
semi-major axis. Reanalysis of the radial velocity data 
indeed reveals such planet closely matching the predicted 
properties. In other words, the eccentricity of the or- 
bit clearly reveals the existence of a perturbing planet, 
and a grazing transit will be most sensitive to even mild 
non-coplanarity between the two objects. We find such 
possible inclination change and we identify a radial ve- 
locity signal matching the predictions of our dynamical 
integrations. 

The resulting planet system around GJ 436 strength- 
ens the trend of a relatively large number of hot- 
Neptunes and supe r-Earths around stars of low mass 
(jBonfils et alJ[20f)l . The presence of a long-term ra- 
dial velocity trend of ~1 m s _1 per year could still be 
indicative of further planets in the system with wider 
orbits. Indeed, the system around GJ 436 shows strik- 
ing resemblances t o that around the M-type star Gl 581 
(|Udrv et alj l2007f ). and thus its planets may experi- 
ence changes in the orbital elements, perhaps eventually 
undergoing transits in spit e of a previously null result 
(|L6pez- Morales et al.ll2006l ). Our study provides yet an- 
other illustration of the variety of exoplanet systems and 
highlights the potential for complex dynamical histories 
that imply sizeable variations of the planets' orbital ele- 
ments over timescales of decades. 

The method of using near-grazing transits should be 
of special interest to discover small planets. Even for 
mild non-coplanarity, objects as small as a few M 9 can 
lead to moderately high values of the induced eccentricity 
and to orbital inclination changes in the order of tenths 
of degrees per year. This will be especially effective for 
transit search missions from space, which could overcome 
the lower detection probability of near-grazing transits, 
and push the detection limits to even lower-mass objects 
than those responsible for the transit events. 
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